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Abstract. We show that a set of conformally invariant equations 
derived from the Fefferman- Graham tensor can be used to construct 
global solutions of vacuum Einstein equations, in all even dimen- 
sions. This gives, in particular, a new, simple proof of Friedrich's 
result on the future hyperboloidal stability of Minkowski space- 
time, and extends its validity to even dimensions. 



1. Introduction 

Consider the class of vacuum solutions to the Einstein equations 
',g) in n + 1 dimensions, which are future asymptotically simple, 
i.e. conformally compact, in the sense of Penrose, to the future of a 
complete Cauchy surface (^,7). A natural method to try to construct 
such space-times is to solve a Cauchy problem for the compactified, 
unphysical space-time (^#,5), and then recover the associated physical 
space-time via a conformal transformation. However, a direct approach 
along these lines leads to severe difficulties, since the conformally trans- 
formed vacuum Einstein equations form, at best, a degenerate system 
of hyperbolic evolution equations, for which it is very difficult to prove 
existence and uniqueness of solutions. 

Friedrich [20, 23] has developed a method to overcome this difficulty 
in 3 + 1 dimensions, by introducing a system of "conformal Einstein 
equations" whose solutions include the vacuum Einstein metrics and 
which transforms naturally under conformal changes. A variation upon 
Friedrich's approach, again in 3 + 1 dimensions, has been presented 
in [6]. 

In this paper, we develop a different approach to this issue which, be- 
sides its simplicity, has the advantage of working in all even dimensions. 
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The method, carried out for vacuum space-times with A > in [2], is 
based on use of the Feffer man- Graham (ambient obstruction) tensor Ti, 
introduced in [16]. The tensor Ti is a symmetric bilinear form, depend- 
ing on a metric g and its derivatives up to order n + 1, cf. Section |2] 
for further discussion. It is conformally covariant, (of weight n — 1) and 
metrics conformal to Einstein metrics satisfy the system 

(1.1) H = 0. 

When n = 3, i.e. in space-time dimension 4, the Fefferman-Graham 
tensor is the well-known Bach tensor. 

The main result of the paper, Theorem 14.11 is the proof of the well- 
posedness of the Cauchy problem for the equation (jl.ljl . for Lorentz 
metrics. This leads to a new proof of Friedrich's result on the future "hy- 
perboloidal" stability of Minkowski space-time [21] (see Theorem 16.1(1 . 
and extends the validity of this result to all even dimensions. 1 As a 
corollary we additionally obtain existence of a large class of non-trivial, 
vacuum, even dimensional space-times which are asymptotically simple 
in the sense of Penrose, see Theorem 16.21 

We further note that in [2] existence of solutions of the Cauchy prob- 
lem for (|1.1|) is obtained by pseudo-differential techniques. Here we 
show that the Anderson-Fefferman-Graham (AFG) equations (|1.1|) can 
be solved using an auxiliary, first order, symmetrisable hyperbolic sys- 
tem of equations. This shows that (|1.1|) is a well posed evolutionary 
system, directly amenable to numerical treatment. Thus, in space-time 
dimension four we provide an alternative to Friedrich's conformal equa- 
tions for the numerical construction of global space-times [18,19,28]. 

Our methods do not apply in odd space-time dimensions, where the 
situation is rather different in any case, as one generically expects poly- 
homogeneous expansions with half- integer powers of 1/r, where r is, 
say, the luminosity distance, compare [11,26,27,30]. 

2. The Anderson-Fefferman-Graham equations 

Let, as before, n + 1 denote space-time dimension, with n odd. The 
Fefferman-Graham tensor Ti is a conformally covariant tensor, built out 

^Once most of the work on this paper was completed we have been informed of 
the work by Schimming [32], who has done a local analysis, related to ours, of the 
Cauchy problem for the Bach equations in dimension four. The application of his 
work to global issues, for instance concerning the vacuum Einstein equations, seems 
not to have been addressed. We are grateful to R. Beig for pointing out that reference 
to us. 
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of the metric g and its derivatives up to order n + 1, of the form 

(2.1) Ti = (V*V)^- 2 [V*V(P) + V 2 (trP)] + T n , 
where 

(2.2) P = Rie 9 -|%, 

and where J- n is a tensor built out of lower order derivatives of the 
metric (see, e.g., [25], where the notation O is used in place of Ti). It 
turns out that T n involves only derivatives of the metric up to order 
n — 1: this is an easy consequence of Equation (2.4) in [25], using the 
fact that odd-power coefficients of the expansion of the metric g x in 
[25, Equation (2.3)] vanish. (For n = 3, 5 this can also be verified by 
inspection of the explicit formulae for and J- 5 given in [25].) 
The system of equations 

(2.3) n = 

will be called the Anderson-Fefferman-Graham (AFG) equations. It has 
the following properties [25]: 

(1) The system (|2.3(l is conformally invariant: if g is a solution, so 
is (f 2 g, for any positive function ip. 

(2) If g is conformal to an Einstein metric, then (|2.3j) holds. 

(3) Ti is trace-free. 

(4) Ti is divergence-free. 

Recall that Ti was originally discovered by Fefferman and Graham [16] 
as an obstruction to the existence of a formal power series expansion for 
conformally compactifiable Einstein metrics, with conformal boundary 
equipped with the conformal equivalence class [g] of g. This geomet- 
ric interpretation is irrelevant from our point of view, as here we are 
interested in (|2.3|) as an equation on its own. 

3. Reduction to a symmetrisable hyperbolic system 

Let g be a Lorentzian metric, let V be a connection (not necessarily g- 
compatible), and let □ denote an operator with principal part g^ u V v 
(acting perhaps on tensors). Let ti be a tensor field, and let V^u denote 
any tensor formed from the A:-th order covariant derivatives of u. For 
k > consider the system of equations 

(3.1) a k+1 u = F(x, u, Vu, V (2) n, . . . , V (2fc+1) u) , 

for some smooth F. Here we allow the coefficients of □ as well as 
the connection coefficients to depend smoothly upon x as well as upon 
the collection of fields («,V«,V^u,...,V' 2f:+1 'u), in particular the 
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(3.3) □ 



metric g is allowed to depend (smoothly) upon those fields. We will 
assume that (|3.1j) is invariant under diffeomorphisms, although this is 
not necessary for some of our results below, such as local existence and 
local uniqueness of solutions. 

We want to show that solutions of 1)3.1 Jl can be found by solving a 
first order symmetric hyperbolic system of PDEs. The idea of the proof 
can be illustrated by the following example. Consider the equation 

(3.2) U 2 u = . 

Introducing 

^(°) = u , = □« , 

it is easily seen that solutions of (|3.2j) are in one-to-one correspondence 
with solutions of the system 

/ v {0) \ _ f 
U (1) J"v o 

It is then standard to write a symmetrisable-hyperbolic first order sys- 
tem so that solutions of (|3.3|) are in one-to-one correspondence with 
solutions of the first order system with appropriate initial data (com- 
pare the calculations in the proof below) . 

Some work is needed when we want to allow lower order derivatives 
as in the right-hand-side of (|3.1|k 

Proposition 3.1. There exists a symmetrisable hyperbolic first order 
system 

(3.4) P$ = H ($) , 

where P is a linear first order operator and H does not involve deriva- 
tives of&, such that every solution of (|3.1|) . with (^,g) time orientable, 
satisfies ()3.4I) , 

Proof. Let {e a } a= o,...,n = {eo, ei}i=i,...,n be an orthonormal frame for g, 
with eo a globally defined unit timelike vector; (such vector fields always 
exist on time orientable manifolds). We set 

(3.5) <pW = Wa}... ai }i<i<2(k~j) , wh ere (p$„. ai = e ai ■ ■ ■ e a p?u , 

(3.6) ip = {(p^} <j< k , 

(3.7) V = {?P {j) }o<j<k , where ^ = CPu . 
Let us derive a convenient system of equations for ip. First, 

(3.8) for 1 < i < 2(k - j) - 1 eo<^U = '^l..a t = M^ 00 ) , 
where we use a generic symbol L to denote a linear map which may 
change from line to line. This gives evolution equations for those 
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i Pai...a i , s which have a number of indices strictly smaller than the max- 
imum number allowed. 

We continue by noting that, again for 1 < i + 2j < 2k — 1, we have 
on the one hand 

n 

(3.9) Q^U = -eo¥'£ 1 ...a i +E e ^£...a,+^ (i) )» 



i=i 



and on the other 
(3.10) Dtp® 



= e ai ---e a n^ + [a,e ai ---e ai }^ 
Combining those two equations we obtain 

n 



Note that the condition i + 2j < 2k — 1 implies j < k so that (|3.11|) can 
also be rewritten as 



(3.12) 



w2...», = E e rf 1 ...«,+ i ^ 



Next, for i + 2j = 2k — 1 and for ^ running from 1 to n we write 
The rewriting of (|3.12j) - ()3.13j) in the form 

+0 



(3.14) 



-eivfi...^ +eorf„ ) ] ... ai +0 



-e \ 



V e nV 3 0o 1 ...a j 



+0 
+0 



()') 

+eo^ai...Oj / 



makes explicit the symmetric character of (|3.12|) - (|3.13|) . It is well known 
that this system is symmetrisable hyperbolic in the sense of [34, Vol- 
ume III] when e° is a nowhere vanishing vector field. 2 This provides the 

(i) 

desired system of evolution equations for those (pa ai—ai s which have 

2 In fact, 113. 1411 is symmetric hyperbolic in a coordinate system with eo = dt and 
eit = 0. However, when g depends upon u and its derivatives it is not useful to 
use such coordinates, as the construction of the Gauss coordinate system leads to 
differentiability loss. In any case Gauss coordinates are not well adapted to the proof 
of existence of solutions when g depends upon u. 
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the maximum number of indices. (One could also use ()3.14j) for any 
number of indices, but 1)3. 8|) is obviously simpler.) 
If we write 1)3.14)1 as 

fV = o, 

where P is a linear first order operator, then the derivatives e a ip satisfy 
a first order symmetrisable-hyperbolic system of equations 

(3.15) Pe a <p = L(ip, Vcp) , L(ip, V<p) := [P, e a ]ip . 
The evolution equations for ip are simply 

(3.16) = , < i < k - 1 , 

(3.17) DV (fc) = F(x,^°\ip(°\Vip^) , 

where in (|3.17|) we have expressed the derivatives of u appearing in 1)3.1 j) 
in terms of (p^ and V<^°) using (|3.5[) . By obvious modifications of the 
calculation starting at 1)3.9)1 and ending at ()3.14)) one can rewrite the left- 
hand-side of ()3.16)) - ()3.17)) as a first order symmetrisable hyperbolic op- 
erator acting on the collection of fields (if>,Vip) := {(V> W , V^ W )}o<i<fe- 
Setting 

(3.18) $ = (^,V^,Vp), 

and letting P be the linear part of the system of equations just described, 
the proposition follows. □ 

The interest of Proposition 13.11 relies in the fact, that it is standard 
to prove existence and uniqueness of solutions of 1)3.4)) when the initial 
data for <J> are in H s , s £ N, for s > n/2 + 1, provided that (^,g) is 
globally hyperbolic. If g does not depend on V 2fc+1 u, then the threshold 
can be lowered 3 to s > n/2. 

Now, not every solution of 1)3.4)1 will be a solution of ()3.1j) . Let 
us show that appropriate initial data for ()3.4j) will provide the desired 
solutions. When the space-time metric g is independent of u, let be 
a spacelike hypersurface in the space-time (.-#, g). We choose eo to be a 
unit time-like vector field normal to «5^, so that the e^'s are tangential at 
J?, and we extend eo off 5? in some convenient way, which might vary 
according to the context. Since 1)3.1)) is an equation of order 2k + 2, the 
associated Cauchy data consist of a set of tensor fields {/(,) }i=o,...,2fc+i 

3 For s > n/2 + 1 the result follows from [34, Volume III, Theorem 2.3, p. 375]. 
However, when the symmetric hyperbolic system has the structure considered here, 
with g not depending upon V 2fc+1 «, the proof in [34] applies for s > n/2. 
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defined on 5? which provide initial data for (ep)*u on S?\ 
(3.19) (eo)^|^ := ^3 u\& = f {i) , < i < 2k + 1 . 

i times 

For any i > and £ > we can use (|3.19|) to calculate formally 

and ip^ \y by replacing each occurrence of (ep^u by f(j\, e.g., |^ = 

f(o), e ^ (0) |j^ = /(1), 

n 

1.5* = (Dtt) I = ( — epep-u. + e^n + r a eo,tt + Tu 



.9> 

n n 

= ~/(2) + e « ei Ao) + r °/(l) + X] ^^/(O) + r /(0) , 
1=1 i=l 

for some linear maps r a , T arising from the detailed structure of □, and 

(i) ■ ( \ (i) 

so on. We will write g)^. for the resulting functions and g)^ for 

the resulting functions epV>^|^, so that 
Similarly we can calculate 



W 1 W 

£ factors 

where we replace each occurrence of (ep) J u by /y), e.g. \y = f^y 

(i) 

We will write h,A for the resulting functions, so that 

V 1^ - >) - % ' V(i)l^ -"(1) ~%) ' ^...Oel^ ~ h (e) ■ 
When g does depend upon u, then the space-time will be built in the 
process of solving the equations. In the simplest case of g depending only 
upon u, the procedure just described should be understood as follows: 
the initial metric g\y is determined by the initial data /p. We choose 
an orthonormal basis {ei}j = o,...,n for g\y, and interpret eo as the unit 
normal to .5? in the space-time that will arise out of the initial data. 
Thus, /n) will be interpreted as the value of the normal derivative of u 
at J?, and so on, and the above considerations remain unchanged when 
this interpretation is used. 

Proposition 3.2. Let $ e C(J, H s {0)), s > n/2 + 2k + 2, s e N, be a 
solution of (j3.4|) on a globally hyperbolic region I x G with initial data 
constructed as described above. Then u := ip^ is a solution of (13. 1|) 
and (EOT?!) . 
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Proof. From (j3.1fij) one has V W = n> (0) for < i < k. It remains to 
show that if (p^... ai = e ai ■ ■ ■ e az ip(°\ then ()3.17|) will coincide with (|3.1j) . 
This can be proved by a standard calculation. One sets 

Xai...a,i = ¥a\...ai ~ e<Jl ' ' ' ea i^^ ' 

and using (|3.4j) one derives a system of equations which show that 

(i) 

Xai-.-oi = for the initial data under consideration. 

However, the computations involved are avoided by the following ar- 
gument. Suppose, first, that g, F and V are analytic functions of all 
their variables. Let us denote by / = {f(j)}o<j<2k+i the initial data for 
(|3.1jl : by an abuse of notation we will write / € H s if £ H s ~ 3 for 
< j < 2k + 1. We note, first, that by using an exhaustion of I X 6 
by compact subsets thereof it suffices to prove the result when I x 6 is 
a conditionally compact subset of the domain of definition of the solu- 
tion. Let f n be any sequence of analytic initial data which converges in 
H s (0") to /. Let u n be the corresponding solution of (|3.4|) ; by stability 
all u n 's will be defined on / x & for n large enough. Similarly, the sta- 
bility estimates 4 for symmetric hyperbolic systems [29] prove that u n is 
Cauchy in C{I,H s (<ff)) n G\l, H 8 ' 1 ^)). The results in [1] show that 
u n is analytic throughout / x & . 

Let u n be a solution of ()3.1|) on an open neighborhood ^ n of in 
1x0" obtained by the Cauchy-Kowalevska theorem. (Note that 9/ n 
could in principle shrink as n tends to infinity, but it is nevertheless 
open and nonempty for each n.) Passing to a subset of ^ n if necessary 
we can without loss of generality assume that ^ ra is globally hyperbolic. 

Now, uniqueness of the solutions of the Cauchy problem for (|3.4|) 
shows that u n coincides with u n on W n . Thus u n satisfies (|3.1|) there 
and thus, by analyticity, everywhere. This shows that maximal globally 
hyperbolic solutions of (|3.1j) with analytic initial data are in one-to- 
one correspondence with maximal globally hyperbolic solutions of (|3.4j) 
with the initial data constructed as above. Then, for H s initial data, 
Proposition 13. 21 follows from continuity of solutions upon initial data for 

Finally, if the fields g, V and F are smooth functions of their argu- 
ments, they can be approximated by a sequence of fields g(n), V(n) and 
F(n) which are analytic in their arguments. The estimates for (|3.4|) just 
described can similarly be used to show that solutions of the approx- 
imate problem converge to solutions of the problem at hand both for 
equation (|3.4|) and 1)3.1 J) , which finishes the proof. □ 

4 Note that I x dff is non-timelike by global hyperbolicity, so that integration by 
parts gives harmless contributions as far as energy estimates are concerned. 
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From what has been said so far we obtain 

Theorem 3.3. Let s > re/2 + 2k + 2, seN. For any fields 

€ H&(S>) , i = 0,...,2fc + l, 

i/iere exists a unique solution of (|3.1|) satisfying (|3.19j) . // i/te metric g 
does not depend upon V 2k+l u, i/ien s > re/2 + 2k + 1 suffices. 

We note that in local coordinate systems (i, x*) on an open neighbor- 
hood % of G C ^ of the form W = I x 0, with ^ n ^ = {t = 0} and 
^-compact, the solutions are in 

tien^Cf/.F-^)) . 

As usual, the lower bound for the local time of existence of the solu- 
tion does not depend upon the differentiability class s, so in particular 
smooth initial data provide smooth solutions. In addition the Cauchy 
problem for (|3. 1|) is well-posed, in that given a pair of initial data fhy 

f"K which are close in H^(^), then the solutions u 1 , u 2 are also close 

When g does not depend upon u, there exists a unique maximal 
globally hyperbolic subset & of , with 5? being Cauchy for 0, on 
which the solution exists. This is proved by the usual methods. 

In the quasi-linear case one also has the existence of a maximal glob- 
ally hyperbolic development of the Cauchy data, giving a space-time 
This follows from the fact that the domains of dependence for 
the system constructed above are determined by the light-cones of the 
metric g, so that a proof along the lines of [5], (compare [7,8]), applies. 

4. The Cauchy problem for the AFG equations 

The Cauchy problem for (|2.3|) has a similar form to that for the 
Einstein equations. Since the system (|2.3fl is of order n + 1, the initial 
data consist of an n-dimensional Riemannian manifold (^,7), re = 
2k + 1 > 3, with n symmetric two-tensors prescribed on 5? . The 
tensor fields represent the i-th time derivative of the metric g in a 
Gauss coordinate system around 5? . Thus, in a neighborhood of 5? , (or 
more precisely a neighborhood of a bounded domain in S fi ), one may 
write 

(4.1) = -dt 2 + 7 (t), 

where j(t) is a curve of metrics on 5? . Setting eo = — Vt, one has 
(4-2) K® = ^4^=0 = \dh[t)\ t=Q . 
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In particular K = is the extrinsic curvature tensor of .5^ in the 
final space-time (^,g). 

The set (7, . . . , K( n ') is not arbitrary, since the equations 

(4.3) W(eo,-) = 0, 

only involve t-derivatives of g up to order re, and so induce (re + 1) equa- 
tions on (7, .K^ 1 ), . . . ,K( n >), Because ()4.3j) is diffeomorphism invariant, 
this is most easily seen in the coordinates (|4.1|) where go a = —So a , so 
that Q4.3J) only involves t-derivatives of g a b, a,b > 1, up to order re. 
The fact that the constraint equations Q4.3|) are preserved under the 
evolution follows in the usual way from the equation STL = 0. 

To describe the system of re + 1 constraint equations Q4.3JI in more 
detail 5 , the Gauss-Codazzi equations for the embedding y C g) 
are: 

(4.4) P 7 - \K\ 2 + H 2 = i? + 2Ric(e ,e ), 



(4.5) SK -dH = Ric(eo,-), 

where H = trK. In addition, in a Gauss coordinate system (t,x l ) near 
y, the Raychaudhuri equation gives 

(4.6) d t H+ \K\ 2 = -Ric(e ,e ). 

We first point out that the curvature scalar R = R g is determined 
directly by the initial data; this is in constrast to the situation with the 
Einstein equations, where R is determined by the evolution equations 
for the metric. Namely, the left side of (|4.4|) is determined by the initial 
data, as is Ric(eo,eo), by (|4.fi|) . Thus R is determined by 7 and 
for % = 1,2. 

To describe the form of the "scalar constraint equation" 7i fiv n^n u = 
0, n = eo, note that from (|2.2|) one has trP = (re — V)R g /2n. Together 
with (gj) and (|Q|> . this leads to 

□Poo - VoVotrP = (□ + V V ) 9 t PT 

(4.7) +i-(D + (l-re)V Vo)P 7 + ... , 

where stands for terms which contain less derivatives of the space- 
time metric. The t-derivatives of the metric of order 4 cancel out, as 

5 An explicit form of l|4.3|l in space dimension 3 can be found in [32]. The param- 
eterisation of the initial data there is rather different from ours. 
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expected, and one easily finds that the equation H fiu n fl n 1 ' = takes the 
form 

(4.8) A— R^ = Pn , 

where, as before, i? 7 is the curvature scalar of 7, with A = D k Dk the 
Laplace operator of 7. Finally, p n is a functional of (7, K^ 1 ' , . . . , K^) 
which does not involve derivatives of initial data of order n+1, while the 
left-hand-side does. This shows in particular that Q4.3J1 is a non-trivial 
restriction on the initial data. 

It also follows from (|4,7|) that p n will contain terms of the form 

(4.9) iA(trK (n - 1} ) - D k D l K%~ 1} , 

with other occurrences of K^ n ~ 1 ^ there, if any, being also linear with 
at most one space-derivative. Thus, the scalar constraint equation 
Ti^n^n 11 = can be viewed either as a non-linear equation of order 
n + 1 which puts restrictions on 7 in terms of the remaining data, or as 
a second order linear PDE for the trace of K^ n ~ l \ 

One can similarly check that the equation Ti^i = takes the form of 
a linear first order PDE for , with principal part 

(4.10) D\K\f - -tTK^ 9ij ) , 

where D denotes the Levi-Civita connection of 7. 

A set (7, . . . , K^) will be called an initial data set if the fields 
(7, K^ l \ . . . , K^ n >) satisfy the constraint equations (|4.3I) . The collec- 
tion of initial data sets is not empty, as every solution of the general 
relativistic constraint equations solves ()4.3|) . 

The equations (|4.3j) are preserved by the following family of transfor- 
mations, related to the conformal invariance of (|2.3j) . Suppose that the 
data set (7, . . . , K^ n ') arises from a space-time (^,g) satisfying 
(tO) . Then (7, . . . , if<») satisfies (JOJ), and if Q is any strictly 
positive function on then the set (7, . . . , K^) obtained on 5? 
from the metric Q 2 g also satisfies (|4.3[) , For example, if we set 

(4.11) u := , := e • • -e (O)|^ , 

j factors 

where eo is, e.g., a geodesic extension of eo off J?, then it holds that 

(4.12) 7 = w 2 7, K {1} = ojK {1) + w (1) 7 . 

Similar but more complicated transformation formulae hold for , i > 
2, see Appendix^] This leads to a family of transformations preserving 
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(|4.3)l which are parameterized by n + 1 functions u, uj^\ j = 1, . . . , n, 
on which are arbitrary except for the requirement that u > 0. 

Recall that under a conformal transformation of the space-time metric 
we have 

4 4 / 4?1 

(4.13) = n— 9il => i? = n-s=i LR - — — n g n 

It follows from this formula that for any given uj = tt\,y > and a/ 1 ) = 
eo(^)|,y one can choose a/ 2 ) so that solves the linear wave equation 

An 

(4.14) -n„n = m 

(n — 1) 

when restricted to 5? . (This equation is globally solvable in globally 
hyperbolic space-times, but this is irrelevant for the current discussion. 
Note that solutions of Q4.14J1 might sometimes develop zeros; these are 
essential in the analysis of the vacuum Einstein equations). Note as 
remarked above that the curvature scalar R is determined by the order 
2 part of the initial data set. Taking further eo-derivatives shows that 
the remaining cjW's may be chosen so that the conformally transformed 
curvature scalar R, together with its normal derivatives up to order n— 2, 
vanish at 5? . 

Now the conformal and diffeomorphism invariance of (|2.3|) requires 
a suitable choice of gauge in order to obtain a well-posed system. As 
in [2], we use constant scalar curvature for the conformal gauge and 
harmonic coordinates for the diffeomorphism gauge; the treatment of 
the conformal gauge is somewhat different here than in [2]. 

Thus, we require first that 

(4.15) R = . 
If (j1~T5|) holds, then (|2~3|) takes the form 

(4.16) (V*V)— Ric= -J* , 

In harmonic coordinates {y a } = {(r, y 1 )} with respect to the conformal 
gauge (|4.15j) . one has 

Ric afe = -]^g^ v d^d v g ab + Q a b{g, dg) , 

where Q is quadratic in g and dg. Applying (V*V)^ n_1 - )//2 to this and 
commuting (V*V)^ n_1 ^ 2 with the dg terms in Q shows that (|4.16jl has 
the form 

n + l ^ 

(4-17) LV 9ap = ~F^ 

where F™g still has the form (|3.1|l and U g = g^d^dy acts on scalars. 
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Finally, as initial data for the (gauge-dependent) variables go a , < 
a < n, we choose 

90a = Soa, on y. 
The data g a b and d T g a i,, a,b > 1, are determined by the initial data 7, 
. The derivatives d T go a are then fixed by the requirement that the 
coordinates {y 13 } are harmonic when restricted to y, i.e. 

(4.18) OyP = d a g a ? + \g aP g^d a g^ u = on <? . 

The higher derivatives d l T g a p, i < n, on y are then determined by 
the given initial data (7, . . . , K^ n >) and by setting to zero the r- 
derivatives of 1)4.18)1 up to order n — 1. 
These choices lead to the following: 

Theorem 4.1. Consider any class [7, K^\. . . , K^)]) satisfying the 
constraint equations (|4.3jl . with 

(7, K^, . . . , kW) g Hf oc (y) x H°~\y) x ... x H^{y) , 

s > n/2 + n+l, s£N, where (y,*y) is a Riemannian metric and where 
the K® 's are symmetric two-covariant tensors; the equivalence class is 
taken with respect to the transformations of the data discussed above. 

Then there exists a unique maximal globally-hyperbolic conformal 
space-time (./#, [g]) satisfying ()2.3)1 . where [•] denotes the conformal 
class, and an embedding 

i-.y^Ji, 

for which 7 is the metric induced on y by g, with given by 1)4. 2|) . 
One can always choose local representatives of [g] by imposing 1)4.15)) . 
Moreover, the Cauchy problem with such initial data is well-posed in 

H? oc (y) x K~\y) x ••• x H£?{y). 

Remark 4.2. The inequality s > n/2 + n suffices for existence of unique 
solutions in local coordinate patches in Theorem 14.11 One expects that 
the use of local foliations with prescribed mean curvature and space- 
harmonic coordinates as in [4] should allow one to lower the threshold 
s > n/2 + n + 1 to s > n/2 + n in this result. 

Remark 4.3. The conformal space-time (^#, [g]) is smooth if the ini- 
tial data are. Similarly, real-analytic initial data lead to real-analytic 
solutions. 

Proof. Given any initial data set [(7, . . . , K^)\ satisfying the con- 
straint equations, by using the functions uj^\ j = 2, . . . , n from (|4.11[) . 
one can adjust the tensor fields K®, i = 2, . . . , n, so that R g , together 
with its transverse derivatives up to order n — 2, vanish on y, (see the 
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discussion following (|4.14jl ). Note that this holds for any 7 £ [7], so 
that ui > and while fixed, are otherwise freely specifiable. Solv- 
ing (|4.17|) with this initial data, as described in Section |3J one obtains 
a collection of space-time coordinate patches with a solution of (|4.17|) 
there when s > n/2 + n. We recall again that the Cauchy problem for 
(|4.17|) is well-posed in the Hf oc spaces above. 

The argument that these gauge choices are preserved, so that 1)4.1 5 J) 
holds and the coordinates y a remain harmonic in these local space-time 
coordinate patches generated by (|4.17|) . is rather similar to the one for 
the Einstein equations [17]; we give details because of some differences 
in the analytical tools used. 

First recall that 1)4. 17|) can be written in the form 

re-l 

(4.19) LV (R lu ,-V li \ v -V v \ ll ) = -J%,, 

where 

A^ := -±LV . 

As g solves 1)4.19)) . its obstruction tensor equals 
(4.20) 

n— 1 f I n — 1 "-3 

= (-i)— jrv (V/A. + v^-— R g g^)-^— v^i^ 

The divergence identity V^W" = gives then an equation involving A 

and R g : 

(4.21) 

— f 1 ] 

□ g 2 < DX U + Vj,(V M A^ — -^Rg) \ = lower order commutator terms . 

Since TL is trace-free, we further have from ()4.20)) 

n-l 1 

(4.22) D 5 2 ( V %--i? g )=0. 

Because R g vanishes to order n — 2 at S*, and Xj vanishes to order n — 1 
at 5? by the discussion following 1)4.18)1 . the initial data for this equation 
vanish. It follows, for instance from the work in Section 13) that 

(4.23) R g = 2V% . 
This can be used to rewrite 1)4.21)1 as 

n + l 

(4.24) U g 2 A„ = lower order commutator terms . 

In this last equation all commutator terms in 1)4.21)1 that involve gradi- 
ents of R g have been replaced by derivatives of A using 1)4.23)1 . 



Since g satisfies the constraint equation (|4.3|) at 1)4.23)1 and 1)4.20)1 
imply that A,- vanishes to order 11 at so that 1)4.24)1 has vanishing 



ASYMPTOTICALLY SIMPLE EVEN-DIMENSIONAL SPACE- TIMES 



15 



initial data. This system has the form considered in Section El so we 
conclude that A vanishes. Hence, the coordinates y a remain harmonic, 
R g = by l)4.23j) . and so H = as well, as desired. 

The usual procedure, as used in the context of the Cauchy problem 
for Einstein's equations, then allows one to patch the solutions together 
provided s > ra/2 + n+ l. An argument as in [5] leads to a unique, 
(up to diffeomorphism) maximal globally hyperbolic manifold [g]), 
with [g] satisfying 1)4. 16|) . with an embedding i : 5? — > with the 
desired initial data on i(y). 

The well-posedness statement follows immediately from the discus- 
sion following Theorem 13.31 □ 

Remark 4.4. In general, there will not be a global smooth gauge for 
(^#, [g]) in which R g = 0. The local coordinate patches where R g = 
need not patch together smoothly, preserving R g = 0, since the initial 
data u > 0, a/ 1 ) on local space-like slices 5? are freely chosen, and so 
not uniquely determined. 

To see this in more detail, consider for example the de Sitter space- 
time ^# = K x S", with metric 

9ds = ~dt 2 + cosh 2 1 0sn(i). 

This is a geodesically complete solution of the Einstein equations with 
Rg = n(n + 1), and so satisfies (|2,3j) . The linear wave equation (|4.14|) 
with initial data $7 = c\, dt^l = C2 on 5? = {t = 0} has a global solution 
on given by 

^ n — 1 . , f 1 , 

it = : — smhr + ci2 / — ——at + di, 

4 ,/ cosh t 

for suitable d\,d,2- One sees that there are no values of c\ , C2 for which 
£1 > everywhere on so that there is no natural global R = gauge 
for (./#, [gds]) with such initial data. 

Consider for example the solution f2 = — sinht, giving an R = 
gauge in the region t < 0, which does not extend to {t = 0}. For t = —e 
small, the induced metric on S^- e = {t = —e} is the round metric 7,5 
on S n of small radius 5 = 5(e). To obtain an R = gauge starting at 
S^s which extends up to and beyond {t = 0}, one must choose ui to be 
a large constant. This causes a discontinity in the choice of gauge for 
the metric, but not in the structure of the conformal class. 

Remark 4.5. One may also construct the maximal solution (^#, [g]) by 
means of local gauges satisfying 

(4.25) R g = c 
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in place of the scalar- flat gauge (|4.15|) . for any cq G M. The proof of this 
is the same as before, noting from the form of (|4.13|) that given uj and 
it/ 1 ) on J^, one can find u/ 2 ) on 5^ such that R = cq on 5? . 

For example, the de Sitter space-time is a geodesically complete so- 
lution in the gauge R g = n(n + 1). However, it is well-known that 
(^#, gas) conformally compactifies to the bounded domain in the Ein- 
stein static cylinder gE = —dT 2 +gs"(\) where T G (— | , ^). The metric 
gE is of course also a solution of (|2.3I) with R gE = n(n — 1), which is 
thus a globally hyperbolic extension of (^,gds), since T € (—00,00); 
it is easily seen that this is the maximal solution. 

Thus, the choice R gE = n{n — 1) gives a global gauge for [gds]- The 
de Sitter metric itself, with gauge R g = n(n + 1), has conformal factor 
Q relative to gE blowing up to 00 as T —* To obtain an extension 
of [gds] past the range (—§,§) of gds requires a rescaling of the large O 
factor to a factor of unit size. 

5. AFG EQUATIONS VS. ElNSTEIN EQUATIONS 

Consider an initial data set (=5^, 7, K) for the vacuum Einstein equa- 
tions in n+1 dimensions, n odd. Thus (J^, 7) is a Riemannian manifold, 
and the pair (7, K) satisfies the vacuum constraint equations with cos- 
mological constant A € M. Using Einstein's equations one can formally 
calculate the derivatives 

(5.1) K® := ^7«|t=o , 1 < i < n 

in a hypothetical Gauss coordinate system near 5? in which the space- 
time metric g takes the form —dt 2 + 7(i), as in (|4.1|) . This gives 

Proposition 5.1. The initial data set for (|2.3|) so obtained solves the 
constraints (|4.3|) . and any suc/t globally hyperbolic solution of (|2.3|1 
given by Theorem \4-l\ is conformally Einstein. 

Proof. Let a be the associated maximal globally hyperbolic solution of 
the vacuum Einstein equations. Then a also solves Q2.3JI . and the result 
follows from the uniqueness part of Theorem 14.11 □ 

A space-time with boundary (^#, g) is said to be a conformal comple- 
tion at infinity of a space-time (^#, g) if the usual definition of Penrose 
is satisfied; that is, there exists a diffeomorphism <I> from jtft to the 
interior of jM and a function f2 : ^# — > M + U {0}, which is a defining 
function for the boundary 

J := &J 7 , 
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such that 

_g = $*(n- 2 g) . 

A conformal completion g) is H s smooth if g G H s (y$), (in suitable 
local coordinates for yt\ 

A set (y, 7, K, oj, u/ 1 ), . . . , u/ n )) is said to be a smooth conformal 
completion at infinity of a general relativistic initial data set (y, 7, K) 
if (y, 7) is a Riemannian manifold with boundary, with y = y U (9^, 
and with 

: y ^> M 

being smooth-up-to-boundary functions such that a; is a defining func- 
tion for the boundary 

y := &y , 

with 

(5.2) 7 = w 2 7 

on =y . Finally K is a smooth-up-to-boundary symmetric tensor field on 
y such that the equation 

(5.3) f = wif+^7 

holds on (compare ()4.12|) ). We further assume that the functions 
a;®, 2 < i < n, are such that the fields IfW, calculated using the if^'s 
as in (|5.1|) . and the functions uW, defined before (|4.12|) . can be extended 
by continuity to smooth tensor fields on y. 

A conformal completion will be said to be H s if 7 G H s {y) and 
e H'-^y), i = l,...,n. 

The above conditions are clearly necessary for the existence of a 
smooth conformal completion d la Penrose of the maximal globally hy- 
perbolic development of (y,~f,K); we will see shortly that they are 
also sufficient. We emphasise, however, that Equations (|5.2|) and (|5.3|) 
alone, together with the requirement of smoothness of 7 and K are not 
sufficient for the existence of such space-time completions. Indeed, it 
follows from the results in [3] that, in dimension 3 + 1, the requirement 
of smoothness up-to-boundary of the K^'s, i > 2, imposes further con- 
straints on 7 and K. It would be of interest to work out the explicit 
form of those last conditions, analogously to [3], in all dimensions. 

We have the following, conformal version of Proposition ^. II it allows 
one to repeat several constructions of Friedrich (see [19] and references 
therein) for vacuum space-times with vanishing cosmological constant 
in all even space-time dimensions: 
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Theorem 5.2. Let (S^,j,K) be a general relativistic vacuum initial 
data set, A = 0, which admits an H s conformal completion at infinity, 
s > n/2 + n + 1, with s S N and n odd. Then there exists an H s 
space-time with boundary g), equal to the conformal completion at 
infinity of the unique maximal development (^,g) 0/(^,7, K), so that 
g = Q~ 2 g on ^ , and 

Proof. The proof is essentially identical to that of Proposition Con- 
sider the initial data (7, . . . , K^) as constructed at the beginning 
of this section. One can conformally transform them to initial data 

R u ■■= ^7fcj|t=o , 1 < i < n 

in a hypothetical Gauss coordinate system near 5? in which the space- 
time metric g takes the form g = —ds 2 + j(s). Here we use a nor- 
malisation as in the proof of Theorem 14.11 or Remark 14.51 requiring the 
vanishing (or constancy) of Rg. Theorem 14.11 provides a solution of 
this Cauchy problem, while the fact that this solution is conformally 
Einstein follows from Proposition 15.11 

The choice of the conformal factor f2 transforming (^#, [g]) to the vac- 
uum Einstein solution g) of course depends on the choice of gauge 
for [g]; we describe here how Q is determined at least in the natural 
settings corresponding to (|4.25|) . Suppose that 7 = lo 2 j is a geodesic 
compactification of (^,7), so that for x near J?, u(x) = dist^J? , x) . 
Such a compactification is uniquely determined by the choice of a bound- 
ary metric on the boundary 5? . Now the value of a/ 1 ) , at the zero level 
set of u is determined by the initial data, (compare Eq. (3.13) of [3] 
in dimension 3 + 1; an obvious modification of that equation holds in 
all dimensions). Changing time-orientation if necessary, one will have 
u/ 1 ) = — 1 at 5? and we extend a/ 1 ) to a neighborhood of 5? in 5? to 
have the same value. These data determine the compactification of the 
initial data set {S fi , 7, K). 

Hence, the proof of Theorem 14.11 gives a unique local solution g of 
()2.3j) in the conformal class [g] , with Rg = cq, for any given cq , satisfying 
the initial conditions. Let <p = S7(™ -1 )/ 2 , and set g = £l~ 2 g = ip~ i ^ n ~ 1 >g. 
Then ip, (and so O), is uniquely determined in a chart for j$ containing 
a portion of by the requirement that <p solves the linear wave equation 
-4rjn<7</2 — Rg<p = 0, with initial data ip = u, d s (f = co^ = —1 on 5? 

near , (where s is the Gaussian coordinate). 

Given such a solution Q, let J* be the connected component of the 
set {£1 = , dil 7^ 0} intersecting 5?. Since, by construction, Ric g = 0, 
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and g is smooth up to J" , standard formulas for the Ricci curvature 
under conformal changes show that 0, has the usual structure on J^, in 
that Vf2 is null, and V 2 ri is pure trace, on . □ 

Remark 5.3. A version of Theorem 5.2, and its proof, also holds for de 
Sitter-type vacuum solutions of the Einstein equations, where A > 0. In 
this case, the completion is at future or past space-like infinity J*^ or 
J^ _ ; the Cauchy data for 7) at consist of the two undetermined 
terms grfy, gr^ in the formal Fefferman-Graham expansion for vacuum 
Einstein solutions with A > 0. This gives an alternate proof of one of 
the results of [2]. 

The case A < leads to initial-boundary value problems. While it is 
clear that a generalization of Friedrich's analysis of this case [22] should 
exist, precise statements require further investigation. 

6. Applications to semi-global and global stability of 
general relativistic initial value problems in all even 
space-time dimensions 

Let {5? ', 70) be the Poincare metric on the (n + l)-dimensional ball, 
with n > 3 and n odd. Setting 

Ko = 7o, 

the set 70, Kq) is an initial data set for the vacuum Einstein equa- 
tions, denoted standard hyperboloidal initial data. The maximal globally 
hyperbolic development (^#,50) of (S?, 70, -Ko) is given by 

(6.1) g = -dr 2 + r 2 70 , 

for t E M + , with = {r = 1}. This space-time is the interior of the 
future light cone about a point in Minkowski space-time (the "Milne 
universe"). With respect to the standard smooth conformal compact- 
ification of Minkowski space-time as a bounded domain in the static 
Einstein cylinder, 

(6.2) g = - dT 2 + dR 2 + sm 2 R g S n-i {1) , 

one has = 0, 1 < i < n; the hypersurfaces {r = const} correspond 
to the level-sets {T = const}, and so are totally geodesic. As noted in 
Remark 14.51 this choice of gauge is global and satisfies R go = n(n — 1). 

As an example of application of the results of the previous section, 
one now easily obtains: 

Theorem 6.1. Let 5? he an n- dimensional open ball, n odd, and 
consider a general relativistic initial data set (S^,j,K) which admits 
an H s conformal completion at infinity, s > n/2 + n + 1, s E N. 
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Then there exists e = e(n) > such that if the associated data 
(y,[{^,K^,...,K^)}) are e-close in H s x . . . x H s ~ n to the data 
[(70, 0, . . . , 0)]) associated to standard hyperboloidal initial data, 
then the maximal globally hyperbolic development (^#,5) 0/(^,7, K) 
is causally geodesically complete to the future. 

The H s conformal compactification {y& , g) is H s close to (*df, go), 
and extends to a larger H s space-time containing a regular future time- 
like infinity l + for (^#,5). 

Proof. The standard space-time (^,go) has a conformal compactifica- 
tion to a bounded domain T> in the static Einstein cylinder (|6.2|) . where 
V corresponds to the range of parameters T + R £ [0, n], T — R € [0, tt], 
R > 0. Future null infinity is given by = {T + R = it}, T G 
(f ,7t), with future time-like infinity t + = {T = ir,R = 0}. The future 
development of (5? , 70, ifo) corresponds to the domain V + = V (l{T > 

-1 

2 

Clearly, the compactification (D,go) extends smoothly to a neighbor- 
hood P of P as a globally hyperbolic solution of (|2.3[) , The Cauchy data 
for such an extension are an extension of the standard Cauchy data 
(5?, [(70, 0, . . . , 0)]) past the boundary 5?. Similarly, the initial data 
(y, [(7,i? (1) , . . . ,i? (n) )]) for (JK,g) extend in H* past a neighbhorhood 
of the boundary and generate a maximal globally hyperbolic space- 
time satisfying ()2.3|) , By the Cauchy stability associated with 
Theorem 15.21 for e small, the solution is close in IP to (V,g~o), 
and in particular is an if s extension of {y$, ', g), where j% = {f2 > 0} in 
jj£ . This shows that, to the future of ,y, (^#, g) has an H s conformal 
completion, which extends in H s to a neighborhood of J? and t + . This 
gives the result. □ 

Note that in dimension 3 + 1 the mere requirement 

l^,K {l \...,K { ^)] € H s x ... x H s ' n , s > n/2 + n + l , 

regardless of any smallness condition, forbids solutions which have loga- 
rithmic terms with small powers of 1/r in poly homogeneous expansions. 
Thus, (similarly to the results of Friedrich), the above theorem applies 
for non-generic initial data sets only. 

Using Corvino-Schoen type constructions together with the above 
stability result, as in [9], one obtains: 

Theorem 6.2. There exists an infinite dimensional space of complete, 
asymptotically simple globally hyperbolic solutions of the Einstein vac- 
uum equations in all even dimensions n + 1 = 2{k + 1), n > 3. Thus, 
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such solutions are geodesically complete both to the future and past, and 
have a smooth conformal completion at infinity. 

Proof. The Corvino-Schoen gluing technique [12, 13] can be used, as 
in [9,10], to construct static, parity-symmetric initial data on R n , for 
any n > 3, which are Schwarzschild with m^O outside a compact set, 
and which are as close to the Minkowskian data as desired. The result- 
ing maximal globally hyperbolic space-time then contains smooth hy- 
perboloids, close to standard hyperbolic initial data, as in Theorem 16. II 
both in the future and in the past. In even space-time dimensions the 
result follows as in the proof of Theorem 16.11 □ 

We note that all the space-times constructed in the proof of Theo- 
rem 16.21 possess a "complete J^" ; this should be understood as com- 
pleteness of generators of J? in the zero-shear gauge, compare [24]. 6 

One expects the above construction to generalise to initial data which 
are stationary, asymptotically flat outside of a spatially compact set 
(rather than exactly Schwarzschild there). This would require proving 
that the resulting space-times have smooth conformal compactifications 
near l°, (in space-time dimension four this follows from [14,15,33]), and 
working out "reference families of metrics" needed for the arguments 
in [10]. Those results are very likely to hold, but need detailed checking; 
note that one step of [33] requires dimension four, and that the family 
of asymptotically flat stationary metrics in higher dimensions might be 
richer than that in dimension 3 + 1, [31]. 

Appendix A. The infinitesimal invariance group of the 

CONSTRAINT EQUATIONS. 

In this appendix we study the Lie algebra associated to the group 
of conformal transformations preserving the constraint equations (|4.3|) . 
This allows one to derive identities which shed light on the structure of 
those equations. 

Suppose that TC = and that the space-time metric g is rescaled by 
Q 2 , where, in a Gauss coordinate system (i,x*), so that 5? = {t = 0}, 
we have 

V 

for some j > 0, and e > small. In the new Gauss coordinates (t, x l ) 
we thus have 

Q 2 (-dt 2 + gijdx l dx j ) = -dt 2 + g kl dx k dx l , 



-'We are grateful to H. Friedrich for useful discussions concerning this point. 
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leading to 

(A.l) (d t t) 2 - g ke d t x k d t x* = Q 2 , 

(A. 2) d t idit = g ke diX k d t x e , 

(A.3) gij = IT 2 (g ke d i x k d j x e - diidji 

By definition of Gauss coordinates we have (i,x l ) = (0{t),x l + 0(t)), 
and also (t, x l ) = (t + 0(e), x' + O(e)). Inserting this in the equations 
above, matching powers in Taylor expansions we find 

i = t+e—^—P +1 +0(e 2 t 2 ^ 3 ) , x* = x i +e-^-^+ 2 +0(e 2 ^ +2 ) . 
U + !)• U + 2). 

At the right-hand-side of (|A.3|) we have variations related to the fact 

that all the quantities there are evaluated at the point x = x+ex (•)+■ • ■; 

to first order, this produces a Lie derivative-type contribution. Next, 

there are variations related to the fact that t = t + e x (•) + ... Each 

term tfKW /i\ in (twice) the Taylor expansion of g~ij at t = gives then 

a contribution to the right-hand-side of (|A.3j) equal to 

From this we can calculate the coefficients of an expansion in powers of 
t of the right-hand-side of (|A.3|) ; inverting those relations, for j = this 
leads to 

, A4 ^ (l-e(i-2)^« + 0(e 2 ), i^2; 

( J "I K«-!^ + 0(e 2 ), i = 2, 

where «£? denotes a Lie derivative; note that ^D^g is twice the Hessian 
of ip. 

For j > we obtain instead 
(A.5) = K® , < i < j , 

(A.6) A"Ci) = + ei>g , 

v ' 2(j + l)!i! 

J 0(^ 2 ), 1<*^2; 
1 -^ D ^g + 0(e 2 ), i = 2. 

The terms proportional to e in those equations describe the desired 
infinitesimal action. 

We can view Tiofi as functions of the metric, its derivatives, and of 

the symmetric derivatives Dr^^^K^J , with i = 1, . . . ,n. As Ti^ is a 
differential operator in the space-time metric of order n + 1, the possibly 
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(i) 

non-trivial contributions arise from those Dr^ ji^K^ for which s + i < 
n + 1. Differentiating the equations = Ho/j, with respect to e, at 
e = t = for j > one obtains 

(A.8) 

n+l-j 

= E — 0jL -^9hlD tl .u^ 

"2 — -^Di^am^ ■ 

At any point x the derivatives D( u-- g a (p(x) = Df£ lM j\(p(x) can be cho- 
sen independently, which leads to various identities. The simplest one 
is obtained for j = n, then the last two lines give a vanishing contribu- 
tion. We parameterise the i£w's by their trace-free part and by trace, 
obtaining 

(A.9) = - dHo>M . , , s>0; 

this gives a check of the general form of (|4.10j) . 

When j = n — 1 the last line in (jA.8j) is zero again, and we obtain 

s =0 o u {£x...£s) K M 



s =o dD tl .. Ia K K M > 



2 V *^Lan. . ?H n ) 



s=0 *1— *s 



' dK ki 9D t K kt 
The vanishing of the coefficients in front of the second derivatives of ^ 
leads to the identity 



(A ' 10) ° dD^trK^- 1 ) 

consistently with (|4.9j) . 
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